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The basic physics of nucleation in solid single-crystal nanoparticles is revealed by a phase-field 
theory that includes surface energy, chemical reactions and coherency strain. In contrast to binary 
fluids, which form arbitrary contact angles at surfaces, complete "wetting" by one phase is favored 
at binary solid surfaces. Nucleation occurs when surface wetting becomes unstable, as the chemical 
energy gain (scaling with area) overcomes the elastic energy penalty (scaling with volume). The 
nucleation barrier thus decreases with the area-to-volume ratio and vanishes below a critical size, 
and nanoparticles tend to transform in order of increasing size, leaving the smallest particles ho- 
mogeneous (in the phase of lowest surface energy). The model is used to simulate phase separation 
in realistic nanoparticle geometries for Li x FeP0 4 , a popular cathode material for Li-ion batteries, 
and collapses disparate experimental data for the nucleation barrier, with no adjustable parameters. 
Beyond energy storage, the theory generally shows how to tailor the elastic and surface properties 
of a solid nanostructure to achieve desired phase behavior. 
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Despite the widespread use of phase-separating 
nanoparticles for catalysis, "smart materials", and en- 
ergy storage [IH3], their complex phase behavior is just 
beginning to be understood @HE]. The most basic open 
question involves nucleation, which is difficult to observe 
experimentally and beyond the reach of ab initio molec- 
ular simulations. (A 20 nm nanoparticle has ~ 50, 000 
atoms.) It is known that surfaces are important [2J, but 
their precise role is unclear. In bimetallic nanoparti- 
cles, where phase transitions are triggered by changes in 
temperature [7J, anomalous melting point depression in 
silica-gold core-shell nanoparticles persists to unexpect- 
edly large particle sizes [5J. Enhanced interdiffusion in 
gold-silver core-shell nanoparticles cannot be explained 
by size-dependent melting point depression [5]. 

For Li-ion battery nanoparticles, where phase transi- 
tions occur by ion intercalation, the literature is full of 
contradictions. The canonical phase-separating cathode 
material is LiFeP0 4 , which exhibits low power in micron- 
sized particles |10j but can achieve very high rates in 
nanoparticles |11) . Experimental measurements of the 
critical overpotential to initiate lithiation vary widely 
from 2mV to 37 mV [T2"U21| . Size dependence has also 
been reported [TJl [21] . Some experiments observe elec- 
trical signatures of nucleation and ID growth above a 
critical particle size [55], while others attribute voltage 
hysteresis to mosaic phase separation among homoge- 
neous particles without nucleation [16 . Some theoret- 
ical studies suggest that nucleation at surfaces leads to 
"intercalation waves" (moving phase boundaries) at low 
current [2"3"H2"5] , while others describe a "solid-solution 
pathway" without the possibility of nucleation [5S] . 

In this article, we resolve these discrepancies by show- 
ing that nucleation in single-crystal nanoparticles is size- 
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FIG. 1: (a) Partial wetting by fluid displacing fluid a 

at a solid surface with contact angle 9. (b) Complete 

"wetting" by solid f3 displacing solid a at a fluid 

surface. Table: Surface energies of FeP0 4 (a) vs. 

LiFeP0 4 (/3) calculated from first principles [55]. The 

smaller phase boundary energy [2T, 7 Q/ g = 39mJ/m 2 , 

implies complete wetting by a or /3 at each facet. 



dependent, occurring as a result of surface adsorption 
that leads to coherency strain (a long-range force). This 
mechanism implies a nucleation barrier that decreases 
linearly with the area-to-volume ratio, in quantitative 
agreement with a wide range of experimental data. We 
show that the data are consistent with both phase-field 
[27] and ab initio calculations |26| that estimate a zero- 
current overpotential of about 35 mV in bulk LiFeP0 4 . 
The analysis presented in Methods is very general and 
could be applied to other multiphase nanostructures. 

Surface wetting in binary solids. - The wetting of a 
solid surface by binary fluids is one of the most stud- 
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(a) Equilibrium voltage curve for a binary intercalation particle, calculated with our phase-field model for 
4 . For Acj) > Acj>* (or A/i > A/i*), transformation from metastable to stable homogeneous states is nucleated 



FIG. 2 
LiFePO_ 

by changes in potential, composition, or temperature leading to unstable two-phase coexistence at the barrier, (b) 

Sketch of nucleation at wetted surfaces due to changes in composition, (c) Calculation of the surface adsorption 

profile in a LiFeP0 4 nanoparticle (C3 shape [2S] with wetted side facets; inset) by voltage fluctuations, which 

increase the bulk concentration, cq(A(J>), until the surface layer becomes unstable at the coherent solubility limit Cq. 



ied problems in fluid mechanics 30-32] . but a theory of 
surface "wetting" in binary solids, which have coherency 
strain, has not been developed. Wagemaker et al. [3"3"ll3~i] 
modeled the effect of surfaces on the voltage curves and 
solubility limits of LiFeP0 4 , but assumed that phase 
boundaries form a contact angle with exterior surfaces 
and neglected coherency strain. Tang and Karma stud- 
ied coherent spinodal decomposition at solid surfaces |35j , 
but did not consider finite systems or dependence of sur- 
face energy on composition. Bai et al. simulated nucle- 
ation by surface wetting [24], but neglected coherency 
strain. 

Young's equation, 7 Q = jp + j a p cos 9, relates surface 
and interfacial tensions to the contact angle 9 at a triple 
junction (Fig. [fja)). Because fluid- fluid and fluid-solid 
interfacial energies have the same order of magnitude, all 
contact angles are possible in binary fluids [22 [22] , in- 
cluding complete wetting and de- wetting (9 = 0°, 180°). 

The surface of a binary intercalation compound in- 
volves equilibrium between two solid phases and a fluid 
(the electrolyte in battery). A stable triple junction 
(0° < 9 < 180°) is unlikely to form between two solid 
phases because the excess energy of the free surface (from 
broken bonds) is much larger than that of a coherent 
solid-solid interface (from stretched bonds) (Fig. fllb)). 
A rule of thumb is that coherent interfaces have j a p < 
200mJ/m 2 , while solid surfaces have 7^,7^ > U/m 2 
[36j . The change in surface energy with composition thus 
dominates, and one solid phase will tend to completely 
wet the surface. If j'(c) = 0, a 90° contact angle will 
form and coherent surface spinodal states exist [35] . but 
this is not the typical situation. 

This prediction is supported by the "core-shell" struc- 
tures commonly observed in bimetallic nanoparticles 
with the low-7 phase as the shell [2J, as well as by 



first principles calculations of battery nanoparticles. The 
change in surface energy by between a — FeP0 4 and 
/3 = LiFeP0 4 for common facets of the Wulff shape 
(Table fl]) greatly exceeds the phase boundary energy, 
lap = 39mJ/m 2 , inferred from the thickness of striped 
phases [27] . As a result, each crystal facet tends to be 
fully lithiated (9 = 0) or delithiated (0 = 180°). 

Nucleation at Solid Surfaces. - Consider a nanoparti- 
cle of volume V in bulk phase a with a surface area A 
wetted by molecules of phase j3, in metastable equilib- 
rium with diffusional chemical potential A/i = up — fj, a , 
where A/i = corresponds to chemical (or mosaic) two- 
phase coexistence. In the case of lithium ion intercalation 
in batteries (6] EH |27] , Acj) = Afx/e = cp e - cp is the in- 
terfacial voltage (of electrons cj> e relative to ions <fi), and 
Acj) = corresponds to mosaic phase separation across 
a set of homogeneous particles [THl |37]. For a neutral 
species, A/i can be controlled by adjusting the reservoir 
concentration of that species. 

The transformation from the metastable phase a to 
the stable phase /3 can be nucleated by fluctuations in 
composition, potential or temperature (Fig. [2|a)). The 
nucleation barrier corresponds to the onset of (unstable) 
two-phase coexistence. At a critical potential, A/i* = 
eA</>*, corresponding to the coherent solubility limit, the 
barrier for composition fluctuations vanishes (AG b = 0), 
and the critical state is an unstable surface layer (Xp — 
0). 

With battery nanoparticles, it is straightforward to 
control the voltage and observe a sudden current as- 
sociated with phase transformation. In this case, the 
bulk solid remains homogeneous at the bulk concentra- 
tion co(A(/>) until the coherent solubility limit is reached 
at the critical potential Acj)* (Fig. J2jc)). The Gibbs free 
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FIG. 3: Phase-field simulation of galvanostatic lithiation of a Lix FePC>4 nanoparticle (C3 shape [29] with a 150 nm 
x 76 nm (010) facet (top)) illustrating nucleation at the coherent solubility limit (X = 0.09), followed by growth 
from the lithiated side facets with interface alignment along (101) planes due to coherency strain. (Movie online) 



energy 



G(c ) = (/(c ) + c Q Ap) Ps V + a(c )A 



(1) 



has contributions from the homogeneous Hclmholtz free 
energy per volume, p s f(co) (where p s is the density of 
intercalation sites), and the excess surface free energy of 
the wetted surface layer, er(c ). Setting G'(cq) — at 
the critical concentration Cq = co(A</>*), we obtain the 
critical potential, 
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The critical potential in an infinite particle, Ap*^ — 
— /'(cq), is the difference between coherent and chemical 
(or mosaic) solubility due to the elastic strain exerted by 
the wetted surface on the bulk solid, which scales with 
volume V, as shown in Methods. Since surface energy 
scales with the wetted area A, the critical potential de- 
creases with particle size, L = V/A, and vanishes be- 
low a critical value, L c = a }f,°L , which corresponds to 

PsJ {Co) 

Ap* = (Fig. pi . Below L c there is no barrier for trans- 
formation. 

Phase-field model. - To more precisely determine Acp^ 
and L c , we use the phase-field method, extended to in- 
clude electrochemistry [6l [24] , coherency strain [27] and 
external surfaces [30]. Phase- field models accurately de- 
scribe both homogeneous and heterogeneous nucleation 
[381 I39j . For an inhomogeneous single-crystal nanopar- 
ticle, the free energy is a functional of the intercalated 
ion concentration c(x, i) and elastic displacement field 
u(x,t), 
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where 7(c) is the surface energy as a function of sur- 
face concentration, k the gradient energy coefficient, 
€ij = 1^- the strain tensor, and Cijki the elastic stiff- 



ness tensor. Overpotential is defined as 77 = 
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the phase-field theory of electrochemical kinetics [5J H7J . 



Dynamical equations are given in Methods, and at equi- 
librium G is minimized with ?/ = 0. 

Here we focus on the effect of surface energy, described 
by the natural boundary condition, 
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which is necessary to enforce SG = [30l [39] . For the 
case of a binary fluid in contact with a solid, Cahn [50] 
showed that Eq. [4]is a generalization of Young's Law. For 
a typical binary solid (Fig. [l| , the left hand side of Eq. [4] 
dominates, leading to complete wetting or dewetting of 
each facet, depending on the sign of 7'(c). This can be 
seen by substituting Eq. [9] and performing the graphical 
analysis of Cahn [3D] . 

The properties of the critical point are derived in Meth- 
ods by defining a as the excess surface free energy of the 
adsorption layer (Fig. [2lc)) and analyzing the equilib- 
rium conditions ^ = and ^- = 0. The derivation 
shows that a is independent of system system size, lead- 
ing to Eq. 2 for the critical potential, jj^- is approxi- 
mated as 
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where Bq is the elastic contribution to interfacial energy 
(scaling with CV,), Cg is the concentration at the coherent 
solubility limit, and X — y J v cdV is the state of charge 
of Li x FcP0 4 . For a very large system, a small amount 
of surface adsorption will have a negligible effect on X. 
As the particle size decreases, however, surface adsorp- 
tion significantly changes X, which in turn creates an 
increased energy contribution from coherency strain. 

Equation [l] is reminiscent of the free energy of forming 
a nucleus in classical nucleation theory (CNT), AG = 
VAG V + -/A, but fundamentally different. CNT is an 
approximation applied to a growing nucleus of a few 
nanometers, while Eq. [2] is rigorous and applies to much 
larger nanoparticles with fixed size. The intercalation 
analogue of the classical critical radius is a binary inter- 
calation particle at the critical size L c where A<fi* = 0. 
Bulk free energy will be dominant in particles larger than 
L c , but the phase state of particle smaller than L c will 
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FIG. 4: Size dependence of nucleation in LiFeP0 4 nanoparticles. (a) Calculated discharge curves of circular 

cylinders of various diameters, (b) Calculated critical potential versus wetted area-to-volume ratio for different 

particle shapes, collapsing onto a master line, (c) Measured critical potentials versus reported particle sizes, falling 

on or below the theoretical limit for a perfect crystal (outside the gray region). 



be controlled by surface energy. Below the critical size, 
the nanoparticle will exist in the homogeneous state that 
minimizes surface energy. 

Since it is based on non-equilibrium thermodynam- 
ics [BJ, our model can be used to predict the dynamics 
of nucleation and growth, driven by surface reactions. 
Phase transformation at constant current with a fixed 
transformation time (C rate) is illustrated in Fig. [3] 
At constant voltage, the transformation time depends 
mainly on the thermodynamic driving force, or overpo- 
tcntial beyond critical voltage, and only weakly on the 
particle size. For small driving force, the particle initially 
lingers near the critical point. The transformation then 
accelerates exponentially due to the release of elastic en- 
ergy (linearly growing overpotential from the tilted volt- 
age plateau), but a detailed theory is beyond the scope 
of this paper. 

Application to LiFeP0 4 . - Parameters for the phase- 
field model were estimated previously by fitting a reg- 
ular solution model for /(c) and a gradient energy to 
experimental solubility data [27] . Based on the fitting, a 
chemical solubility limit of X = .01 and coherent solubil- 
ity limit of X — .09 were estimated, in agreement with 
experiment. Using these phase-field parameters and the 
surface energies in Table [I] (Eq. Q reduces to Dirichlet 
boundaries with c s = or c s = 1 with the application 
of Eq. [9]), we consider realistic LiFeP0 4 particle geome- 
tries characterized by Smith et al. |29) . Simulation of 
intercalation dynamics is presented Fig. [3] for particle 
shape C3, using the 2D depth-averaged model described 
in Methods [BJ [57J . The lithiated phase originates from 
the side facets, and intercalation waves [23] HI] propa- 
gate inward, with phase boundary orientation tending to 
occur along elastically preferred (101) planes. 

Figure |4[a) shows calculated equilibrium discharge 
curve for LiFeP0 4 nanoparticles with circular cross- 
sections of varying diameter. The critical potential on 
discharging, or lithiation from a low state of charge, is 



reduced in smaller particles, as a result of lithium adsorp- 
tion on the side facets, which competes with coherency 
strain to facilitate nucleation. (The same effect is not 
seen during charging, or delithiation from high concentra- 
tion, since the 2D depth-averaged model does not model 
dewetting of the (010) surface.) 

The size effect resulting from surface adsorption is cal- 
culated for a variety of particle shapes in Fig. Bb) . Each 
LiFeP0 4 geometry is treated as a prism, with the lithi- 
ated side facets contributing to the wetted surface area 
A. The calculated critical potentials are close to a mas- 
ter line given by Eq. [21 so the A/V ratio is much more 
important than the precise shape. As shown in Methods, 
each surface layer is localized and planar, so its tension 
(energy/area) is largely unaffected by the bulk geometry. 

At a critical wetted area-to-volume ratio, A/V = 
L^ 1 s=s .18 /nm, the critical potential vanishes, A<fi* = . 
Below the critical size L c (analogous to the critical radius 
in classical nucleation theory), the nanoparticle becomes 
dominated by its surface properties, rather than its bulk 
properties. For prism-shaped LiFeP0 4 particles, this cor- 
responds to mean particle diameter, d c = AL C — 22 nm. 
A delithiated particle below this size will spontaneously 
lithiate and remain lithiated at equilibrium. Interest- 
ingly, this is the same size at which phase separation is 
suppressed in LiFeP0 4 [27J HO]. It is also roughly the 
cutoff size for anomalous phase behavior in bimetallic 
nanoparticles [8] [9] . 

The calculations in Fig. Efb) are compared with ex- 
tensive experimental data in Fig. Hfc) without adjusting 
any parameters. For each experiment, A/V is found by 
approximating the particles as prisms and taking A to 
be the area of the side facets. Since the calculation rep- 
resents an ideal thermodynamic limit, experimental data 
points are expected to lie on or below the theoretical line 
(Eq. §. 

The collapse of experimental data in Fig. W[c) strongly 
supports the theory. No experimental points lie signifi- 



cantly above the simulated limit, and most lie very close 
to it. Some points also lie below, as expected for two 
reasons. Firstly, it is likely that some experiments over- 
estimate the size of the active particles. If the system 
has a distribution of particle sizes, a consequence of Eq. 
[2] is that small particles will be transformed first at lower 
potentials. Therefore measuring an average particle size 
may not be sufficient to characterize a porous electrode. 
If measurements are made while only a fraction of the 
total capacity is cycled (as done by Dreyer et al. 16 ), 
the properties of the smallest particles in the system may 
inadvertently be measured. Secondly, defects may play 
an important role by reducing coherency strain and de- 
creasing the barrier for intercalation. Meethong et al. 
measured different amount of retained strain in differ- 
ent samples [TH [T3]. If the amount of retained strain 
decreases, Eq. [2] predicts the bulk phase transformation 
barrier A<C will also decrease. In the limit of zero strain, 
the particle will transform at the miscibility gap where 
A^ - 0. 

Discussion - Li x FeP0 4 was initially thought to be a 
slow-rate material due to kinetic limitations of nucle- 
ation and growth [10]. However, this hypothesis is in- 
consistent with the rapid rate improvements from smaller 
nanoparticles [3] , doping to improve electrical conductiv- 
ity [JTJ , and the use of exotic surface coatings [321 [13] • 
These modifications are not expected to significantly al- 
ter growth kinetics. 

Recently there has been a shift in thinking to the oppo- 
site extreme, with ultra-fast rate capabilities attributed 
to particles that never phase-separate. Malik and Ceder 
[26] calculated a solid solution pathway using quantum 
Monte Carlo and argued against the possibility of nu- 
cleation and growth by applying CNT at the bulk spin- 
odal. A significant result of phase-field theory, however, 
is that the nucleation barrier goes to zero - and CNT fails 
catastrophically - at the spinodal point, in agreement 
with experiment |38j . The argument was also based on a 
phase boundary energy of 7^3 = 960 mJ/m 2 from ab ini- 
tio calculations j33| , which is likely more than an order of 
magnitude too large. We have inferred j a! 3 = 39 mJ/m 2 
from experimental morphology and solubility data [2"7] . 
consistent with the bound 7 Q ^ < 200 mJ/m 2 for coher- 
ent interfaces [36] . (From a modeling perspective, this 
illustrates the limitations of ab initio calculations that 
are restricted to a few unit cells of material, often at 
zero temperature, although our results also demonstrate 
the remarkable accuracy of the ab initio bulk and surface 
properties, used to parameterize our phase-field theory.) 

Dreyer et al. [IS] also advocated for particles that 
never phase separate based on the observation of a 20 mV 
(A</>* = 10 mV) "thermodynamic" hysteresis Li x FeP0 4 . 
Although the basic picture of mosaic phase separation 
may hold, we have shown that there is no unique value 
of the critical potential. In order to be consistent 
with experimental observations of coherent phase sep- 



aration [37], the bulk critical potential must be much 
larger than 10 mV, and the inferred value Ai^ = 35 mV 
is confirmed here by an independent analysis of nucle- 
ation (Fig. El). Indeed, the discrepancy over hysteresis 
can be resolved by combining our theory of nucleation 
with porous electrode theory based on non-equilibrium 
thermodynamics [37] . as will be reported elsewhere. 

Conclusion - We have developed a quantitative phase- 
field theory of nucleation in single-crystal nanoparticles. 
A key observation is that complete "wetting" by one solid 
phase is typically favored at each surface, so that nucle- 
ation corresponds to the instability of a wetted surface 
layer. In order to overcome coherency strain, the nucle- 
ation barrier becomes a linear function of the surface- 
to- volume ratio, which implies that nanoparticles tend 
to transform in order of increasing size. The theory is 
confirmed by collapsing disparate experimental data for 
nucleation in Li x FeP0 4 nanoparticles, without any ad- 
justable parameters. Beyond important applications to 
Li-ion batteries, however, the theory also has broader 
relevance for nanotechnology. It provides the basic prin- 
ciples to design solid nanostructures with desired phase 
behavior under different operating conditions, by control- 
ling the elastic and surface properties of the component 
materials. 



METHODS 

Dynamical Model. - The mean intercalation rate at the 
surface is related to the local surface concentration, stress 



state, and overpotential 77 



SG 



by the phase-field theory of 



electrochemical kinetics [6]. At equilibrium, G is minimized, 
and 77 = 0. To model intercalation dynamics (Fig. |3J), we 
assume reaction limitation for anisotropic LiFeP0 4 nanopar- 
ticles (with fast diffusion and no phase separation along the 
[010] crystal axis and negligible diffusion along the [100] and 
[001] axes) and solve the electrochemical Allen-Cahn reaction 
equation P l23l l2il |27] . 



dc 2 a / 5G 

— - = -J (c, V c,ii)smh — 

at e \ep s oc 



+ i 



(6) 



for the depth-averaged concentration c(x, y) over the active 
(010) facet. The right side of Eq. p]is a generalized Butler- 
Volmer rate for symmetric electron transfer, where Jo is the 
exchange current per area [6]. Stochastic intercalation is mod- 
eled with Langevin noise £, which facilitates nucleation or 
spinodal decomposition [241 i27j. The strain field is deter- 
mined by mechanical equilibrium, ji = V ■ a = 0, and zero 
surface traction, n ■ a = 0, for a solid particle in a liquid 
electrolyte [251157] . 

Derivation of the Critical Potential. - For functionals that 
do not depend explicitly on x, the Beltrami Identity is an inte- 
grated form the the Euler equation that applies at equilibrium 
in ID systems. Application of this relation to Eq. [3]produces: 



p s [/(c) + eA^c] 



-k(Vc) + -CijkiUjtki — C (7) 



where C is a constant whose value can be determined by con- 
sidering the boundary condition at x = L, which is stress-free 



and c(L) — cq. Thus C — p a [/(co) + eA<f>co], where Co is the 
bulk concentration in equilibrium with Acj>. Solving Eq. [7] for 
the gradient energy produces: 



-k(Vc) = p s [A/ + eA^Ac] + -Cijkidjeki 



(8) 



where A/ = /(c) — /(co), and Ac : 
Vc produces: 



Co- 



Solving Eq. [8jfor 



Vc: 
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The excess free energy of the adsorption layer a 30 (which 
differs from the surface energy 7) is the difference in energy 
between a homogeneous system and a system with the ad- 
sorption layer: 



a = G - f Ps [/(co) + eA0c o ] dV 
Jv 



(10) 



Substitution of Eq. 15] to eliminate the gradient energy leads 
to: 

a = 7 (c s ) + 2 / p s [A/ + eA(j>Ac] + -C ijk i6ij6 k i dx (11) 
Jo * 

which is minimized at equilibrium. Eq. [9] can be used to 
change the variable of integration in Eq. |ll|from x to c: 



a — 7(c s ) +-/k \Jlps [A/ + eA^Ac] + Cijkitijtu dc 



(12) 
For an adsorption layer at a surface in a semi-infinite sys- 
tem, the elastic energy can be approximated as Cijkifyjtkl = 
Bq{c— Co) 2 |44| . The stability of the adsorption layer is related 
to the sign of ^p, and a is maximized when ff = 0: 






2([p s /'(c)+eA0] +Bo(c-c ))=0 (13) 



The equilibrium condition eA</> = — /'(co) can be applied in 
the bulk, far from the interface: 



Psf'(c) + B c = p s f'(c ) + BqCq 



(14) 



which is the common tangent construction for coherent binary 
solids H^j. Thus the adsorption layer will become unstable 
and grow at the coherent solubility limit (miscibility gap). 
The energy needed to reach the coherent solubility limit is 
the barrier energy to initiate transformation of the particle. 

The dependence of the surface excess a on the size of the 
system can be determined by examining the sign of |^ , which 
is found using the fundamental theorem of calculus: 

-^ = 2 P s [A/ + eA<f>Ac] + C i]k ie XJ e k i\ (15) 

OLi \x—L 

Since the boundary condition at x = L is stress-free and 
c(L) = c : 

Thus a is independent of the system size L. 
From Eq. |10[ we can see that: 



G= [ Ps [/(c„) + eA0c o ] dV + oA 
Jv 



(17) 



The terms inside the integral are constants, and the integral 
can be evaluated: 

G = Ps [/(co) + eA<M)] V + a A (18) 

Now set 4^- = for Co = Cq and solve for Acj>*\ 

^= P s{f'(c Q ]) + eA<l ) ]V+^(co)A = (19) 

dco aco 



eA<j>* 



-'™ + 7.(&*>)v (20) 



where Cq is the concentration at the coherent solubility limit, 
determined by the coherent common tangent construction 
(Eq. [Ml. 

Using Eq. |12| an approximation for elastic energy 
Cij k i€ij€ k i — Bo(c — c) 2 [44], and applying the fundamental 
theory of calculus, an estimate for ^- can be obtained: 



da 

dc,; 



(c ( *) = sfnB {c* - c) 



(21) 



where c = i J c dV is the mean value of c. 
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